Exact solutions of Brans-Dicke cosmology with decaying vacuum density 



A. E. Montenegro Jr.^ and S. Carneiro-'^'^ 
^Instituto de Fisica, Universidade Federal da Bahia, 4^0210-3^0, Salvador, BA, Brazil 
'^International Centre for Theoretical Physics, Trieste, Italy* 

We investigate cosmological solutions of Brans-Dicke theory with both the vacuum energy density 
and the gravitational constant decaying linearly with the Hubble parameter. A particular class of 
them, with constant deceleration factor, sheds light on the cosmological constant problems, leading 
to a presently small vacuum term, and to a constant ratio between the vacuum and matter energy 
densities. By fixing the only free parameter of these solutions, we obtain cosmological parameters in 
accordance with observations of both the relative matter density and the universe age. In addition, 
we have three other solutions, with Brans-Dicke parameter ua = —1 and negative cosmological term, 
two of them with a future singularity of big-rip type. Although interesting from the theoretical 
point of view, two of them are not in agreement with the observed universe. The third one leads, in 
the limit of large times, to a constant relative matter density, being also a possible solution to the 
cosmic coincidence problem. 



I. INTRODUCTION 

Recent observations suggest that the total energy den- 
sity in the universe is greater than the (barionic -I- dark) 
matter density. On one hand, dynamical estimations 
lead to a ratio between the matter and critical densities 
around one third On the other hand, measurements 
of anisotropies in the microwave background radiation 
indicate that our universe is spatially flat 2], suggesting 
the existence of another, unknown form of energy, usually 
called dark energy. Its presence is also corroborated by 
age estimations [3] and by the observed distance-redshift 
relation for supernovas la, which suggests that, in the 
present phase of universe evolution, the deceleration pa- 
rameter q is negative, and therefore the universe performs 
an accelerated expansion 

Several models have been proposed in order to explain 
those observational data. There are some candidates for 
dark energy, as, for example, the cosmological constant, 
the so called quintessence, or the generalized Chaplygin 
gas. Among them, the simplest and oldest one is the cos- 
mological constant, also associated to the energy density 
of vacuum. 

We can contextualize the study of dark energy in dif- 
ferent gravitational theories. For instance, the theory of 
General Relativity, where the Einstein field equations are 
used, and where the gravitation constant G is a universal 
constant. Another one is the Brans-Dicke theory 0,0], a 
scalar-tensor theory in which the gravitational constant is 
a function of space-time, and where a new parameter, w, 
is introduced. More recently the interest on this kind of 
theory was renewed, owing to its association with super- 
strings theories, extra-dimensional theories and models 
with inflation or accelerated expansion 0]- [Hi- 
lt is generally assumed that General Relativity is re- 
covered in the limit lo —>■ +oo (despite the existence of 
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Brans-Dicke solutions for which this is not true [12] )■ As- 
tronomical observations in the realm of Solar System im- 
pose a very high inferior limit for lo. Nevertheless, such 
a result corresponds to the weak held limit, and applies 
only in the simplest case of constant lo. Therefore, it is 
possible that General Relativity is not adequate to de- 
scribe the universe at early times, or needs corrections in 
the cosmological limit. 

In this article we consider the Brans-Dicke theory and 
associate to dark energy the equation of state of vacuum. 
We investigate models in which the vacuum energy den- 
sity decreases with the universe expansion, a hypothesis 
that has been considered as a possible solution to the 
cosmological constant problem, that is, to the question 
of why the presently observed value of A is about 120 or- 
ders of magnitude below the value predicted by quantum 
field theories 

Our goal is to find solutions of Brans-Dicke theory 
which satisfy a particular variation law for G. We shall 
use the Eddington-Dirac relation, based on the large 
number coincidence, G « H/m'^, where H = d/a is 
the Hubble parameter and m-jr is the pion mass [iH, [l^ . 
We will then take G = H/8ttX, where A has the order 
of m^. In addition, as usual, we will relate the Brans- 
Dicke scalar field to the gravitational constant through 
(j) — Ga/G, where Go is a positive constant of the order 
of unity. 

Together with that variation law for G, we shall con- 
sider two different ansatzen. The first one is given by 
p = 3aH'^ /8ttG, where p — pm + Ph is the total den- 
sity, and a is an adimensional constant of the order of 
unity. This ansatz is suggested by observations, which 
show that Pm ~ Pcl'i, where pc = iH'^ /8tiG is the criti- 
cal density. On the other hand, we know that pA has, at 
most, the same order of magnitude as pm, otherwise its 
presence would be more evident. This ansatz was already 
considered in [13, Ull • 

The second ansatz will be given by A = (3H^ , where (3 
is a constant of the order of unity. We are, in this case, 
inferring a variation law for the cosmological term, which 
has already been considered in the literature on the basis 
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of different arguments [Ta|-[23| (for other variation laws 
for the cosmological term, see, for instance, [23]-[13)- We 
will show that this ansatz leads to a set of solutions larger 
than the first one, containing its solutions as a particular 
case. 

We will look for solutions for recent times, that is, 
we shall consider a spatially flat (fc = 0) Friedmann- 
Robertson- Walker space-time, filled with a perfect fluid 
whose matter component is pressureless (pm — 0). For 
the cosmological term, we will take the equation of state 
of vacuum, pA = — Pa- 



II. SOLUTIONS WITH VARYING 
COSMOLOGICAL TERM 



A. The first ansatz 



Taking pm — 0, k — and pA = — PA, the Brans-Dicke 
equations Q can be written as 



dt 



3 + 2w 



(/9 + 3pA)a^, 



Sirp 



7H 



"6 02- 



(1) 



(2) 



(3) 



We then have a system of three ordinary differen- 
tial equations, with four unknown functions of time: a, 
Prm PA and (p. The system becomes solvable if we add 
the Eddington-Dirac relation, G = H/8ttX, which re- 
lates a and cj). In order to restrict our class of solu- 
tions, we will take in addition our first ansatz, given by 
p — SaH"^ /8ttG. In this way, we obtain 



3aXH, 



and 



= 87rAGo(l + g), 



(4) 



(5) 



(6) 



where q = —aa/o? is the deceleration factor. 

With the help of equations ([H)-®, we can rewrite ([T|)- 
([3]) in the form 

(3 -f 2oj)\GQ[q + (1 + q)iH] = 3aXH + 3pA, (7) 



Pm = aXH{l + q), 



^^2 + ,--(l + ,f. 



(8) 



(9) 



Equation ([9]) tell us that q is constant, since a, Gq and 
uj also are. Therefore, q = 0, and equation ([7|) reduces to 



(3 -I- 2w)AGo(l + q)H = aXH + pA- 



(10) 



By using ^ and ([5]), we can obtain the vacuum den- 
sity. 



PA = aXH(2 - q). 
Substituting A = H/SttG, we obtain 

a{2 - q)H^ 



(11) 



PA 



SttG 



Since pA — A/SttG, we conclude that A — a{2 — q)H^, 
which suggests our second ansatz, with = q;(2 — g), to 
be used later. 

Leading (jlj and fS]) into ((2]), one obtains 



1 

H 



(l + q)t + G, 



where G is an integration constant. Let us take C = 0, 
in such a way that iJ — > 00 for t — + 0. We then have 



H 



1 1 



(12) 



l + qt 

Substituting d/a for H in ([72]) . we also have 

a = At'^, (13) 

where A is another integration constant. 

The relative density of matter, defined with respect to 
the critical density, can be obtained by using G = H/8ttX. 
Then, p^ = 3XH, and, using ([5]), one obtains 



Pm 
Pc 



a(l + q) 



(14) 



Substituting pA from (fTT|) into ^TU\i . we have as well 
{3 + 2Lu){l + q) 



a 
Go 



3-9 



(15) 



Comparing a/Go given by equations (O and p3|) . one 
can derive a relation between lo and q, given by 



-il + qf\{3-q). (16) 



(3 + 2w)(l + g) 



Eliminating uj from equations ^ and (jl5p . we can also 
obtain a/Go as a function of q only, 



_ 12{2 + q) + 3{1 + q)^ 
Go ~ {1 + q){3 - q) + 12 ■ 



(17) 



With these relations, it is easy to derive some results 
to be compared with current observations. For example, 
if g = 0, from equations (fT6|) and (fT7|) we obtain u; — 6/5 
and a/Go = 9/5- From we have Ht — 1. From 
equation (fT5|) it follows that a = At. And, from (|14p . 
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one has f2m/a = 1/3. Since a « 1, we can see that 
^lrn ~ 1/3, in agreement with astronomical estimations 
The age parameter Ht is also in good accordance 
with globular clusters observations [3]. 

If, on the other hand, we would take g = — 1, we would 
obtain, instead of equation (fT^ . the result H ~ constant, 
that is, the de Sitter universe, with /9,„ — and a constant 
PA. Note, however, that q — —1 does not satisfy equation 
p6)) . that is, the de Sitter universe is not solution of 
Brans-Dicke equations for this ansatz. 



B. The second ansatz 

Taking now our second ansatz, A = /3H^ and G — 
H/8ttX, and reminding that pA = A/SttG, we obtain 

PA = I3XH. (18) 

Furthermore, as well as in the first ansatz, we have 

SttXGq 



and 



H 



(/) = 87rAG'o(l + g). 



(19) 



(20) 



With the help of IHIj-JlOl), we can put (P)-® in the 
form 

(3 + 2cu)XGo[q + 3H{1 + g)] = p + SPXH, (21) 



p + 3Hp - 3I3XH^ = 0, 



p = 3XGoH 



2 + g--(l + (?)2 
6 



(22) 



(23) 



Here we have a solvable system, with three differential 
equations for three unknown functions, H, p and q. By 
finding p one can, using (HH]), determine pA and pm- 
Leading p given by ((23|) into pTjl . we derive 

£_ _ (3 + 2uj)[q + 3H{1 + q)] - 3H[2 + g - |(1 + gf] 
Go 3H 

(24) 

In this way, since f3/Go is constant, there are two pos- 
sibilities: either g = 0, in which case we have a sim- 
ple relation between g and /3/Go, or p4)) is an evolution 
equation, with g ^ 0. 



1. The case g = 



In this case, equation becomes 

^^{3 + 2uj){l + g)-\2 + g~^{l + qf 
Go Lb. 



(25) 



Using ([23)) into ([22)) . one obtains, by integration, 
1 3Go[2 + q-^il + qf]-B 



H 



Go[2 + g-|(l + g)2 



Let us choose C — 0, such that H 
we have 



where we have defined 

Go[2 + g-|(l + g)2] 



-t-hG. 



oo for t ^ 0. So 



3Go[2 + g-|(l + g): 
Substituting d/a for in equation 
a = At"\ 



we find 



(26) 



(27) 



(28) 



where A is an integration constant. On the other hand, 
from (pS)) we can obtain q = (1 — n)/n, or 

1 



(29) 

i- -r g 

Leading n from equation ([ST)) into ([29)) . one obtains 

(30) 



Go 



2 + g- g(l + g)'J (2-g), 



and we can write 
PA = f3XH ^ XGaH 



6 



(l + qf (2-g). (31) 



By using and (I?!]) , we also have 



p - PA = AGoi? 



9-7(1 + 9)' (1 + 9)- (32) 
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We can also obtain the relative density of matter from 
equations ([50)1 and ([5^ . Reminding that pc = 3A-ff, we 
have 



nr. 



Prr. 
Pc 



P f l + q 
3 \2-q 



(33) 



Eliminating uj from ([25)) and ([50)) . one can obtain /3/Go 
as a function of g only. 



J3_ _ 12(2 + q) + 3{1 + q)"^ 
Go " (l + g)(3-g) + 12 



(2-g). 



(34) 



On the other hand, comparing /3/Go given by ([25)) with 
that given by (|30p . we obtain the same relation between 
UJ and g we have obtained with the first ansatz, equation 

m- 

This is not a mere coincidence. If we compare p given 
by ([^ with that obtained from equation ^ of the first 
ansatz, we obtain the equation ([9]) of the first ansatz. 
Then, ([?!)) can be reduced to equation (jlip of the first 
ansatz. Equation ([5^ . on the other hand, is reduced to 
equation ([5]) of the first ansatz. 
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Equation ([30)1 can be put in the form 



(35) 



already anticipated in the first ansatz. By using it, it is 
possible to verify that equations p3|) and (|34p are the 
same as and ([T7)) of the first ansatz. Finally, one 
can also verify, with the help of (|^^ , that equations 
and are identical to equations (|12p and of the 
first ansatz, respectively. 

We thus conclude that, in the case of a constant q, the 
two ansatzen are equivalent. 



2. The case q / 

In the differential equation ([H]), substituting d/a for 
H and separating the variables, we obtain 



dq 



da 



3 + 2uj a 



Y Go(6^+6)^+6^(Go+/3) ^ _ , 6+Juj\'^' 

We will initially analyze the case in which the quantity 
2 _ Go{6uj + 6)^ + 6uj{Gq + (3) 

is positive. In sections 2.3 and 2.4 we shall analyze the 
cases in which it is negative or zero, respectively. 

With > 0, let us integrate the above equation by 
doing 



Then, we obtain 
6 



z = q + 



dz 



6 + 7w 



Its solution is given by 



A 



K + q 



6+7uj 



da 
a 



K ~ q 



where A is an integration constant. 

By defining B = k+{6 + 7w)/w, C k - (6 + 7uj)/lj 
and D — [3 + 2lu)/ujk, we have 



a^ A 



B + q 



C~q 



Then, 



while 



^ <Q ^ a ^ A' 



C-q 



q-C 



B + q r. JB + q 



C~q 



C-q 



(36) 



(37) 



(38) 



Let us solve equations ((37)) and ((38)) . in order to find the 
functions q, H and p. 



2.1 - The solution of equation 



Introducing x = a/ A and inverting equation p7|) . we 
obtain 



q 



B + CxT^ 



XD — 1 

With the definition of g, it becomes 

x^{B + CxT^) + xx[x'^ - 1) = 0. 
By taking y = x and y' — dy/dx, one has 

y{B + Cx^) + yx{xT^ - 1) = 0, 
which solution is 



(39) 



where Ci is an integration constant. 

It is easy to see that H ~ y/x, leading to 



Hix) 



Cix 



(X* - l)D(B+C) 

On the other hand, we have dt — dx/y, and so 

{XT^ - 1)^(-B+C) 



t 



Cix^ 



dx. 



(40) 



(41) 



Equations ()39p - (|lT)) are solutions of the Brans-Dicke 
equations (|^T)) and (|^5)) . Let us now verify whether they 
satisfy the third Brans-Dickc equation, (|22)) . Introducing 
p' = dp/dx, and using p given by ()23p . equation 
becomes 



3AG'o 



dH r 

dx 



2 + q--il + q)^ 



3XGqH 



dq 
dx 



3-u{l + q) 



- (sXGaH \2 + q-'^il + q)^] - (3Xh} = 
X I I 6 J J 



By using (|5g)l and (|^ . we obtain 

(3 _ {G + l)xi + {B - 1) 
Gq xi — 1 



{B + C)xo 
D{xT> - 1)2 



;(l + g) 



2 + g-i^(l + g)^ 



-(1 + 9)' 
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Substituting in this equation x given by (|37p . we then 
have 

. ^ BC{lj - 3) 



+ 



{B - C){u) - ?,) - BCuj uj' 
iD{B + C) 2 



Lj(B-C)+CJ-3 

iDiB + C) 



- 1 



In this equation, w, /3/Go, -B; 



315(5 + C) 6 

■(42) 

C and D are constants. 



Therefore, for a varying q, the coefficients of g and 
must be identically zero, simultaneously. This is only 

possible for oj ~ —1 and /3/Go = —3. As Gq is positive, 

we conclude that /3 is negative, that is, the cosmological 

term is negative. 

We then have B = Vl2 + 1, G = \/T2 - I and L» = 

— In this way, equations (l39l) - (|4T|) may be written 

as 



VT2 + 1 + (\/T2- l)a;-^ 



r-VT2 _ I 



H = Cix 



Cit = 



VT2 



(43) 



(44) 



(45) 



In the last one, we have taken a second integration con- 
stant in such a way that a ^ when t — > 0. 

Solution pS]) is plotted in Figure 1. From it we note 
that equation (|37p originates two different universes. In 
one of them, we have a = at t = and, when t — > +oo, 
a/ A 1 asymptotically. On the other hand, q — \/T2 — 1 
for a = 0, tending to +oo as t — s- +oo. Since q is positive, 
the expansion is decelerated, with its velocity tending to 
zero as a/A ^ 1. 

In the second universe, we have the origin of time in 
— oo, with a expanding from its asymptotic value a/ A = 1 
to +00, when Cit — l/\/l2. On the other hand, q is 
initially — oo, when a/ A —t 1, and increases approach- 
ing — (vT2-|- 1) asymptotically, as a/ A -\-oo. As q is 
negative, the expansion is always accelerated. 



0.5 



-2 



FIG. 1: Solution (gSJ; a/ A versus Cit. 



H = 



Cix 



B~l 



D(B+C) ' 



(47) 



and 



t = 



D(B+C) 



■ dx. 



(48) 



with B, C and D defined as before. 

Equations (gGll-dlll) are solutions of ^ and Let 
us verify whether they satisfy the third field equation, 
([22]) . Once more, performing the same transformations 
used to solve (|57)) . we arrive at the same equation ([^^ 
derived before. As we have seen, only with w = — 1 
and /3/Go = —3 we have the coefficients of q, q^ and 
q'^ identically zero, simultaneously. Then, we have again 
B = \/T2 + 1, G = \/T2 - 1, and D = -1/VT2. 

Therefore, equations ((38|) and (|46p - (|48|) can be written 

as 



q = 



Vl2 + l + q 
\/T2- 1 - q 



iVl2- l)x-^-VT2~ 1 



/12 



H = C,x^ (x-^ + 1 



(49) 



(50) 



(51) 



2.2 - The solution of equation 



By performing the same transformations and the same 
steps used to solve (|57|. we obtain the following equa- 
tions: 



Cx^ -B 

XD + I 



(46) 



Cit 



^12 (a 



VT2 . 



1 



12 



(52) 



In the last one, we have chosen the second integration 
constant in such a way that a — > when t 0. 

Solution ((52|) is plotted in Figure 2. We have a — 
at i = 0, and a +oo when Cit l/\/T2. On the 
other hand, q — \/T2 — 1 for a = 0, decreases with the 



6 
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0.05 0.1 0.15 0.2 0.25 
FIG. 2: Solution fSl}: a/A versus Cit. 
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FIG. 3: Solution Sl^/Go versus a/A. 



expansion, becomes negative, and tends to — (vl2 + 1 
when a — s- +oo. 

With the help of we can express equations ^51 
and (|5^ as functions of g: 



H = — 



48Ci 



g2 + 2g- II' 



712-1 - q 
24 ' 



(53) 



(54) 



with q in the interval (— \/l2— I,\/T2— !]• 

The age parameter, on the other hand, can be obtained 
from ([55]) and ([M]) . leading to 



Ht 



q + ^/Vl + l' 
By using HH]), ^ and ([531), we obtain 



Pn 



-24AGoCi 



+ 8g + 19 
g2 + 2(J- II 



(55) 



(56) 



The relative density of matter can then be obtained with 
the help of equations (|53|) and ((56)l . and is given by 



19). 



(57) 



Equations (|56|) and ((57|) can be expressed in terms of 
x, by using ([50|) . We have 



(4 - Vl2)a;^ 



12 



(4- 



I2)a;' 



/12 



and 



Go 



(4 + V12)2; 



12 



(a;-v^ + I 
which is plotted in Figure 3. 



(58) 



(59) 



2.3 - The case < Q 



The solutions found in the previous sections were de- 
rived from the differential equation by assuming 
> 0. Let us now suppose that it is negative. We 
can solve equation ([M]) by doing 



1^1 



Go(6^ + 6)2 + 6u;(Go + /3) 
Goc^2 



and 



6 + 7w 



In this way, we obtain 



6 



dz 



2, + 2uj 



UJ J Kf + Z-' 

which solution is 



n I 6 + 4c^ q 
a — Li exp arctan — 

UJKl 



da 
a 

6+7lj 



where Gi is an integration constant. 

Introducing Di = ujki/ (6 + 4a;) and E — {6 + 7uj)/uj, 
one has 



a = Ci exp arctan 

Dl Kl 



(60) 



Taking its inverse function and defining x = ln(a/Gi), 
we obtain 

q — —Kitaii{Dix) — E, (61) 
or, by using the definition of q, 

X + x^[\ - Kl ta,n{Dix) ~ E]^Q. 
After doing y — x and y' = dy/dx, we obtain 
y' + y[l - Kl tan(£)ix) - £;] = 0, 
whose solution is 

y{x) = = G2e(^"i'^[cos(Dia;)]~"ST, (62) 
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where C2 is a second integration constant. 
As dt — dx/y, we also have 



later is satisfied only if g = 0, in which case equation ([65l) 
reduces to (PO]) . already studied. 



[co,spix)]^i 



(63) 



2.^ - The case 







Equations ([6T|l -(|63 p are solutions of Brans-Dicke equa- 
tions (PTjl and As before, let us verify whether they 
also satisfy the remaining equation, (|22[) . Introducing 
p' = dp/dx and using p given by equation (|^ be- 
comes 



SACq 



diJ r 

dx 



+ 3\GoH 



dq 

dx 



3-Uj{l + q) 



S^SXGqH 1 + 9-^(1 + 9)2 -/3XH^=0. (64) 



By using ([S^ . one obtains 

= [£: + 2 + Ki tan(i:>ix)] [2 + g ^^(1 + g)^ 
Go L 6 



dq 
dx 



3-Lu{l + q) 



Now, with the help of we arrive at 



since dq/ dx — q/ H . 

On the other hand, from one can obtain 

g _ 2 + g-|(l+g)^-(3 + 2c.)(l + g) + A 
3i7 3 -I- 2w 

Leading this expression into (j65p . we have 



(65) 



3.^ 3C.2 



It is possible to verify that (j42|l and ([66|) are identical. 
In the later, w and /3/Go are constants. Therefore, for 
a varying g, the coefficients of g, g^ and g'^ must be si- 
multaneously zero, which is only possible if = — 1 and 
/3/Go = —3. But, in this case, k{ = —12, contrary to our 
initial supposition that k\ > 0. 

Therefore, in the case g 7^ 0, equation ([50]) satisfies the 
Brans-Dicke equations ^ and (USD, but not (HH). The 



In order to fulfill all the possible cases (and solutions) 
let us suppose that = 0. From we have 



6 f dz 



3-t-2w 



da 



where we have done, as before, 

6 + 7tJ 



Its solution is 



1 



LO a 
■ in ■ 



2(3 -H2w) Gi^ 



where Gi is an integration constant. 

Taking D2 = cj/[2(3 2w)] and, as before, E ^ [& + 
7uj)luj, one obtains 

q + E Gi 
Introducing the new variable x — ln(a/Gi), we then have 

1 



D2X 



E. 



(67) 



With the definition of g, this becomes 

x^ + 1)2(1 - E)xx'^ + D2xi = 0, 
or, by doing F = i:>2(l - E), 

x^ + Fxx^ + D2XX = 0. 
Taking now y — x and y' — dy / dx, we obtain 
y + Fxy + D2xy' = 0, 

which solution is 

Fx -|- In a; ^ 



y{x) — H — C2 exp 



D2 



(68) 



where G2 is another integration constant. 

As in the previous case, equations ((67|) and ((68)) are so- 
lutions of (|2ip and , but we should also verify whether 
they satisfy (|22|) . Introducing p' — dp/dx and using p 
given by equation reduces, as we have seen, 

to (|64|) . Now, using equations ([68)l and ([67| leads to the 
same equation (|65p of the previous case, which, as already 
seen, leads to ()66p . Therefore, solutions with varying g 
are only possible if a; = — 1 and P/Gq = —3. With these 
values, however, = 12, contrary to our initial suppo- 
sition that — Q. Wc thus conclude, also in this case, 
that ([22|l is not satisfied, that is, there is no solution with 
varying g. 
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III. CONCLUSIONS 

In this work we have found some exact solutions of 
Brans-Dicke cosmology, by using two different ansatzen. 
We have shown that the first ansatz is a particular case 
of the second one, when the deceleration parameter q is 
constant. 

In the first ansatz, the ratio between the energy den- 
sities of matter and vacuum is constant, characterizing 
a possible solution for the cosmic coincidence problem, 
that is, the approximated coincidence presently observed 
between pm and pc- This possibility survives to a quan- 
titative analysis, since a relative matter density around 
1/3, as indicated by observations, leads to an age param- 
eter Ht 1, corresponding to a universe age around 14 
Gyr, also in accordance with observational limits. 

Nevertheless, this ansatz presents some problems as 
well. The most severe of them is the presence of a con- 
stant deceleration factor (equals to zero if Ht = 1). In 
spite of the claim of some authors (see, for example, 
[28l . [29j ) about the possibility of a uniform expansion 
along the whole universe evolution, a decelerated phase 
is usually considered necessary for large structure forma- 
tion. For this reason, we should consider this ansatz valid 
only in the limit of late times, restricting in this way the 
predictive power of the model. 

With the second ansatz, besides the case of constant q, 
we have found three other universes, with varying q, in 
which the dark energy density is negative and the Brans- 
Dicke parameter is w = — 1. In one of them the decel- 
eration parameter is always highly positive. In a second 
one, it is always highly negative. Therefore, these two 
cases are interesting just from a theoretical viewpoint. 

In the third case, on the other hand, the deceleration 
parameter is initially positive, becoming negative at later 
times, but always finite. In this case (as well as in one 
of the previous cases) one has a future big-rip, with the 
scale factor, the matter density and the Hubble parame- 
ter diverging in a finite time, but with the relative matter 
density remaining finite. As one can see from equation 
([55)) . for an age parameter in the interval 0.8 < Ht < 1.3, 
as defined by the observational limits, the deceleration 
parameter is —2.0 > q > —2.9. Whence, with the help 
of (|57p. it is possible to see that we have, for the rela- 
tive matter density, 1.1 > ^Im/Go > 0.7. As we know, 
different observations restrict the matter density param- 
eter to 0.2 < il„i < 0.4. Therefore, this solution sat- 
isfies such observations, provided Go is in the interval 
0.3 < Go < 0.4. Furthermore, for the whole evolution we 
have fim/Go < 7.5, tending, for future times, to a con- 
stant value around 0.5 (see Figure 3). This also charac- 
terizes a possible explanation for the cosmic coincidence. 

It is interesting to observe that, in the three cases with 
varying q, the cosmological time varies linearly with q, 



which may, therefore, be used to define the time mea- 
surement. It is interesting to note as well that the total 
energy density can be negative, since the dark energy 
density is negative. In our last solution, for example, p 
is positive until q ~ —2.3, becoming negative since then. 

There is a particular result which may seem a limi- 
tation of our solutions, namely the typical values found 
for the Brans-Dicke parameter w. As we know, obser- 
vations in the realm of Solar System impose very high 
inferior limits for it. Let us remember, however, that we 
are considering the simplest version of scalar-tensor the- 
ories, which plays just an effective role here. Corrections 
to General Relativity, if exist, may be scale-dependent, 
and, therefore, observations in the Solar System cannot, 
in principle, impose limits to corrections at the cosmolog- 
ical scale. Particularly, we should not expect any time de- 
pendence of the Brans-Dicke scalar field (and so of G) in 
the Solar System, where the metric is stationary. While 
no spatial dependence should exist in large scale, because 
of the cosmological principle. 

Anyway, a generalization of the solutions studied here 
seems to be necessary, either by modifying our ansatzen 
in the case of early times, either by considering more 
general scalar-tensor theories, with uj depending on the 
scale. The study performed here, though limited in its 
scope, shows the variability of solutions in these contexts. 

We should also note that, for the solution of Brans- 
Dicke equations, it is enough to add, for instance, the 
Eddington-Dirac relation. The inclusion of additional 
constraints, even when empirically or theoretically justi- 
fied, had the goal of limiting the set of solutions. There- 
fore, a possible line of investigation would be to relax 
our ansatzen, imposing to the Brans-Dicke equations, 
for example, just the constraint given by the Eddington- 
Dirac relation, enlarging in this way the class of possible 
solutions. Such a generalization may also include the 
radiation phase, although we do not have any indica- 
tion about the validity of the Eddington-Dirac relation 
at early times. 

To conclude, let us remind that, despite the analysis 
we have made about the observational limits for the age 
and matter density parameters, a more detailed analy- 
sis of the whole set of current observational data is still 
in order. In particular, a careful study of the distancc- 
redshift relation for supernovas la constitutes the subject 
of a forthcoming publication. 
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